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In this paper, we show that if passive fluctuations are considered, primordial black holes (PBHs) 
can be easily produced in the framework of single-field, slow-roll inflation models. The formation of 
PBHs is due to the blue spectrum of passive fluctuations and an enhancement of the spectral range 
which exits horizon near the end of inflation. Therefore the PBHs are light with masses < 10 15 g 
depending on the number of e-folds when the scale of our observable universe leaves horizon. These 
PBHs are likely to have evaporated and cannot be a candidate for dark matter but they may still 
. . . affect the early universe. 
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Inflation [l| is becoming a standard model for the very early universe. The inflationary scenario, in which the 
present universe is only a small local patch of a causally connected region at early times which underwent an expo- 
nential expansion driven by the inflaton potential, is generally accepted for explaining the observed spatially flat and 
homogeneous universe. In addition, its quantum fluctuations during inflation give rise to primordial Gaussian matter 
density fluctuations with a nearly scale-invariant power spectrum, which is consistent with recent astrophysical and 
Oh cosmological observations such as structure formation and cosmic microwave background anisotropics [2[ . 

' Although the simplest single-field, slow-roll inflation model works well, some basic questions have yet to be an- 
qj . swered. What is the origin of the inflaton potential? Do classical matter density imhomogenities that we observe 

■ today genuinely come from quantum fluctuations of the inflaton? Are the observed matter density fluctuations truly 
Gaussian? How robust are the predictions for a subdominant contribution of tensor modes to the metric fluctuations, 
a slightly broken scale invariance, and a negligible running spectral index of the power spectrum? Future cosmic mi- 
crowave background measurements and mega-scale mappings of the large scale structure will definitely answer some 

■ of these questions or perhaps pose a challenge to the standard inflation scenario. 

There has been a lot of studies on inflationary models that go beyond the simplest single-field, slow-roll inflation. A 
, class of models has considered a new source for generating inflaton fluctuations (so-called passive density fluctuations) 
during inflation through a direct or gravitational coupling between the inflaton and other quantum fields. This leads to 
\Q , very interesting results such as the so-called warm inflation Q , the suppression of large-scale density fluctuations Q , 

■ possible constraints on the duration of inflationary expansion Q. the bursts of particle production that result in 
infra-red cascading Q, the trapped inflation in which the inflaton rolls slowly down a steep potential by dumping 
its kinetic energy into light particles at the trapping points along the inflaton trajectory and electromagnetic 
dissipation in natural inflation [9H12I]. 

In all of these papers, essentially, the generation of passive density fluctuations is originated from quantum fluctu- 
ations in the back reaction of the couplings to the inflaton perturbation. The nature of passive fluctuations is usually 
non-Gaussian and non-scale-invariant. A particular feature is that the power spectrum of the passive fluctuations 
can be very blue l8HT2j . These passive fluctuations cannot dominate the primordial density perturbation at large 
scales as confirmed by cosmological observations such as cosmic microwave background (CMB) anisotropics and the 
formation of large scale structures. However, depending on individual models, their significant contribution to the 
non-Gaussianity is still possible. This will be tested soon in the Planck CMB mission and in future large-scale- 
structure surveys. In this paper, we point out that the passive fluctuations with a blue spectrum can dominate the 
primordial density perturbation in the very small scales, seeding the formation of primordial black holes (PBHs) in 
the radiation-domination era after inflation. 
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II. PASSIVE DENSITY FLUCTUATIONS DURING INFLATION 

Let us consider a slow-rolling inflaton <f> coupled to a certain quantum field x- The Lagrangian that is relevant to 
is given by 



L=-g^d^d u( j ) -V{4>)+L I ^,x), (1) 



where V is the inflaton potential, Li is the interaction term, and the metric is 

ds 2 = dt 2 -a 2 (t)dx. 2 . (2) 
The equations of motion for the mean fields arc then given by 

2 B _ , 



+ 3^ + y' = ^, (4) 

where p x is the energy density of x particles, the dot and the prime denote differentiating with respect to time and 
4> respectively, and Mp = 2.4 x 10 18 GcV is the reduced Planck mass. In Eq. (j4]), the righthand side of the equation 
is the back reaction of the interaction to the inflaton mean field. The back reaction arises due to copious production 
of x quanta during inflation. 

In addition, the fluctuations of <f> satisfy 

V 2 f 81 \ 

5<f> + 3HS<f>-^S<l> + V , '5<t> = 8\-^J. (5) 

The homogeneous solution of this fluctuation equation gives rise to the standard primordial density fluctuations. Here 
we call them as active fluctuations and denote their power spectrum by P a . The righthand side of Eq. ([5]), which comes 
from the fluctuations of the backreaction, acts as a source for generating additional fluctuations of </>. The particular 
solution of Eq. ([5]) with the source term is referred as passive fluctuations and the power spectrum is denoted by P p . 
Hence, the total power spectrum is given by the contributions from both active fluctuations and passive fluctuations 
as 

P = Pa+Pp= (S-) (\SM 2 ), (6) 



where 6<pk is the Fourier mode of 



III. PRIMORDIAL BLACK HOLES 



For single-field slow-roll inflation, the spectrum from active fluctuations is given by 



where the slow-roll parameter is 



P = 1 V - (7) 

a 24^ 2 A/4 e' () 



Mp (V 



In a given region with radius r, the criteria of black hole formation is given by 

2GSM 2SM 



rM 2 



>1 (9) 
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where G is Newton's constant and 5M is the mass inside the region r. 1 The condition can be expressed by using the 
energy density 5M ~ Spr 3 as 



M 



Sp>^- (10) 



The primordial density perturbation can be imagined as density fluctuation between different Hubble patches of the 
universe (so called 'separate universes' (HI) with radius r ~ 1/H. Namely, in each patch of the universe, the energy 
density is regarded as homogeneous but each patch has different value of energy density. Therefore Eq. (|10[) becomes 
Sp > AIpH 2 ~ p where in the second equality Friedmann equation (p ~ 3H 2 Mp) is used and p should be regarded 
as the average of many patches. From this simple estimation, we can naively guess that if P 1 / 2 ~ (Sp/p) > Oil), the 
whole patch of the separate universe will colapse into a black hole. The argument here is heuristic one. For those 
who concern about the gauge dependence of density fluctuation, our argument has assumed a spatially flat gauge in 
which we choose a slice with zero curvature perturbation and hence the relevant quantity is density fluctuation. More 
rigorous calculation for the condition of primordial black hole formation is given by Refs. [l4|, HH as 

\<s= 6 i<i- (ID 

3 p 

The upper bound is to avoid formation of a separate closed universe. If the spectrum is P 1 ' 2 ~ 0(1O~ 2 ), PBHs will be 
copiously produced [Til ]. This is certainly not the range of the matter power spectrum probed by CMB experiments. 
The range corresponds to the scale of quantum fluctuations that exits the horizon during inflation at a number of 
c-folds, N = 60, before inflation ends. It is because we have P(N = 60) = (5 x 10~ 5 ) 2 from CMB observations. 
However, the spectrum can be large near the end of inflation, 

P(N = 0) - 10~ 4 - 10 6 x P(N = 60). (12) 

In this case black holes will form soon after inflation when the scale enters the horizon. These black holes are called 
primordial black holes (PBHs) (see Refs. [I(l[l2| for review). The spectrum from active fluctuations is (almost) scale 
invariant, therefore unless the running spectral index is large, PBHs cannot be formed (l8l. [l9|. We can see from 
Eqs. (O and (J5J that the slow-roll parameter e has to decrease toward the end of inflation in order to enhance the 
(active) spectrum. However, inflation has to end so the usual tendency is an increasing e. This is the reason why it 
is so difficult to have PBHs (for single- field slow- roll inflation with active fluctuations) j20l - |23j . However, we point 
out for the first time that it is very natural and easy to have PBHs formed even in the case of single-field slow-roll 
inflation if we consider passive fluctuations. 

IV. INFLATION MODELS AND PASSIVE POWER SPECTRA 

In this section, we will discuss the passive fluctuations in different inflation models and explore the possibility of 
forming primordial black holes from the passive fluctuations. 

A. Axion Inflation 

Let us consider the case that the inflaton is a pscudo Nambu-Goldstone boson with a typical potential after the 
shift symmetry is broken, 

V{4>) =A 4 [l-cos(^//)], (13) 
where A is a mass scale and / is the axion decay constant. And the coupling to a gauge field is 

Li = —^F^. (14) 



1 Interestingly this relation can be found from Newtonian physics by requiring the escape velocity to be larger than the speed of light. 
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This coupling is natural in the sense that there is no symmetry to forbid it. It is found [1(1 [H| ( and been further 
analysised in Ref. [l2[ ) that the inverse decay of the gauge field will enhance the fluctuation of the inflaton field. The 
spectrum is given by 

P = P a fl + 7.5x yr*pj^\ , (15) 



£ 6 , 

where the second term corresponds to the passive spectrum P p and 

^2H] = y2—- (16) 

At CMB scale (TV = 60), large £ would produce large non-Gaussainity according to Eq. (fl~5]) . Therefore, there is 
an upper bound from CMB data which is roughly £ < 3 [h], HH- However, £ becomes larger after CMB scale exits 
horizon. For example, if we approximate Eq. (|13|) by a quadratic potential V = m 2 ip 2 /2 with m = A 2 // which is 
a good approximation when cf> <C /, 2 we would have the slow-roll parameter e ~ 0.01 at N = 60 (CMB scale) and 
e ~ 1 at N = 0. This means that £ grows ten times larger according to Eq. (|T6|) . The total spectrum (near the end 
of inflation) is plotted in Fig. Q]by using Eq. (|T5|) . We can see from the figure that when £(N = 0) ~ 6, Eq. (fTTj) 
is satisfied and PBHs would form. This would correspond to £(N = 60) ~ 0.6 (noting that this can be achieved by 
a ~ 1 if / < Mp), which is well below the constraint (£ < 3) from Refs. pil hj. This means that we do not have 
large non-Gaussianity but we have PBHs. Actually the abundance of produced PBHs may put an upper bound on 
the non-Gaussianity generated in this model. 

One may concern about the effect of backreaction if we have £ ~ 6 near the end of inflation. First of all, this effect is 
milder in our case than the case where large non-Gaussianity is generated at CMB scale. In both cases, £ is expected 
to grow toward the end of inflation and we consider a smaller value of £ at N ~ 60. It has been shown in Ref. [l2j 
that the backreaction to the homogeneous inflaton field due to particle production slows down the inflaton velocity 
and prolongs the end of inflation (instead of shortening the period of inflation). However, this does not change our 
results since the power spectrum depends only on £, which in turn depends on the slow-roll parameter e, and the end 
of inflation is determined by e = 1. 

B. Inflation with Trapping 

Consider an interaction between <p and massless scalar \i °f the- type, 

L I = c l 2 Y J {^~^) 2 xl (17) 

i 

where g is a coupling constant and fa is a constant field value. When <f> rolls down to each trapping point at 4>i along 
the inflaton trajectory, the Xi particles become instantaneously massless and are produced with a number density that 
increases with </j's velocity. As <f> dumps its kinetic energy into the Xi particles, it is slowed down and the produced 
Xi particles are diluted due to the inflationary expansion. 

It is useful to define a time scale, At = l/yjgcj). Then, We have 

HAt= (- P a \ (18) 

For HAt < 1 or g 2 > 10" 7 , it was shown [y, LZ| that bursts of Xi particle production takes place in a time scale, At, 
and the number density of the Xi particles produced is given by n Xi ~ 1/At 3 . This leads to the backreaction of the 
Xi particle production to the inflaton field in Eq. @ and the righthand side of Eq. © is particle number density 
fluctuations in the process of particle production. Furthermore, for closely spaced trapping points with equal spacing 
r along the inflaton trajectory, trapped inflation occurs even on a potential which is too steep for slow-roll inflation, 
provided that HT/(j) < 1 0- 



2 Actually this will result in / > Mp which is not phenomcnologically desirable, but it can be evaded for example if we consider the 
potential as an effective potential of N-flation [24j . See Q| for more examples in multi-field or extra-dimensional models where it is very 
natural to have / < Mp and a/f S> 1/Mp. 
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In the case for HAt ~ 1 (i.e. g 2 ~ 10" 
Extrapolating the result obtained in Ref. | 
given by 



and r ~ </>Ai, it was found [8[ that P p has a blue power spectrum, 
we find that the total power spectrum in this case is approximately 



P = Pa 



1 + 0. 



"CMS 



(19) 



where the pivotal scale kcMB = 0.002Mpc~ , if inflation lasts for about 60 e-folds. This reproduces the power 
spectrum in Figure 2 of Ref. [|[ for k/kcMB < 300. Therefore, P can easily reach unity for small scales that exit 
horizon before the end of inflation. Here we have assumed that the scale dependence does not change dramatically 
during inflation when doing the extrapolation. This is plausible since inflation is (quasi) de Sitter phase and hence 
no extra physical effect is expected to change this scaling behavior during inflation. This blue-tilted behaviour is 
reasonable because a certain fc-modc is being kicked constantly by the trapping fields all along the inflaton trajectory 
and the kicks result in accumulative growth of the fluctuations (see Eq. (54) of Ref. Q ) . The growth is bigger for 
higher k modes since they leave the horizon later in time. 



C. Quantum Stress Tensor Fluctuations 



The last model does not involve a direct non-gravitational interaction between inflaton and other quantum fields 
as discussed above. The passive fluctuations come from gravitational effects and are therefore intrinsic and model 
independent. In Ref. the authors considered the effect of quantum stress tensor fluctuations of conformal fields 
such as conformally coupled masslcss scalars and electromagnetic fields, denoted by ST^, in dc Sitter spacetime upon 
the expansion of a congruence of timclikc geodesies. They found that the quantum stress tensor fluctuations induce 
a source term on the righthand side of Eq. (|5|) , given by 



(20) 



V 

— 69, 

where 9 is the expansion of the congruence, related to the covariant derivative of the fluid four-velocity by 
In the presence of quantum stress tensor fluctuations, they found that 

59=- a - 2 (t)[ dt'a 2 {t')5R )llv ui*uv, (21) 

J t 

where t is when inflation begins and 5R^ V is the Ricci tensor fluctuations induced by 5T^ V , from which the passive 
power spectrum was obtained as 



P P = 10 



-40 



S 



E 



10 12 GcV 



k 



3Mpc" 



(22) 



where S is the total expansion factor during inflation and Er ~ V 1 / 4 is the reheating energy scale at the end of 
inflation. For example, if we consider the scale leaving horizon near the end of inflation at which PBHs form, we 
would have k ~ e 60 kcMB- Assuming that S > e 60 and Er = 10 12 GeV, we can easily obtain a condition that P p ~ 1. 



PRIMORDIAL BLACK HOLES AS DARK MATTER? 



When a PBH forms, roughly speaking the whole horizon collapses into a black hole and therefore the mass is given 
by the horizon mass, 



M 



4tt 



= 47T- 



H 



0.066M, 



© 



— Y 

GeVy 



(23) 



where the last equality assumes radiation domination (with energy density pn = (ir 2 /30)gT 4 ) when black holes form 3 , 
Mq ~ 10 33 g is the solar mass, T is the temperature, and g is the statistical degrees of freedom. By using the formula 



3 For estimation and simplicity, we assume that reheating happens immediately after inflation and therefore PBHs always form in a 
radiation-dominated universe. 
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FIG. 1: Power spectrum produced near the end of inflation. If £ ~ 6, primordial black hole is formed. 

of entropy density, s ~ gT 3 , the conservation of the total entropy, S = a 3 s, implies that T oc g~ x l 3 ja. Thus, the PBH 
mass can also be written as a function of comoving wave number k as [25| 



M ~ 6.4 x 10 14 M 



(-) 

\50J 



-1/6 



Mpc~ 



(24) 



For example, if the PBH forms at T = 10 8 GeV, we have M ~ 10 15 g. A PBH with mass smaller than 10 15 g would 
have evaporated through Hawking radiation [26| . PBHs with larger masses contribute to (or might even dominate) 
the dark matter density. The mass fraction /3 = Pbh/p of PBHs of mass M is given by [13] 



/3(M) 



dS 



1/3 V2ttct{M) 



exp 



2ct 2 (M) 



<j{M) exp - 



18cr 2 (Af) 



(25) 



where er(M) is the mass variance at the horizon crossing. The range of integration corresponds to Eq. (fTTj) . If PBHs 
form during the radiation-dominated epoch, a ~ 0.655. Hence, the density parameter Q BH (M) of the PBHs in the 
present universe is [25j 



n BH (M)h 2 ~ 2.1 x 10 8 /?(M) 



(Gev) ' 



(26) 



where /i is the present Hubble constant in units of 100 km/scc/Mpc. For example, if the PBHs with M ~ M Q play 
the role of dark matter in the present universe, i.e. ilsHh 2 ~ 0.25, then /? ~ 5 x 10 -9 , a ~ 0.06, and 5 ~ 0.092. 
For the scale of PBH formation leaving horizon near the end of inflation, we would have k ~ e 60 kcMB ~ 10 23 Mpc — . 
By using Eq. (|24p. we have the PBH mass M ~ 10 _32 A/© ~ 20g. These small black holes would have already 
evaporated and could not be dark matter. We thus conclude that PBHs generated from scales leaving horizon near 
the end of inflation are unlikely to become dark matter. However, if we consider low-scale inflation in which the CMB 
scale corresponds to N ~ 45, we may have PBHs with M ~ 10 15 g that are evaporating now. This may explain the 
observed antiproton fluxes from the BESS experiment [ID, Furthermore, PBHs with M < 10 15 g can have other 
cosmological implications. For examples, they could affect baryogenesis (29j and nucleosynthesis |30L l3l| . swallow 
monopoles [H, [3^], destroy domain walls and so on (see Ref. [35| for more detailed references). 



VI. CONCLUSION 



A lot of effort has been put in studying the passive fluctuations during inflation and their imprints on the CMB and 
in particular the non-Gaussian features. However, in this paper we have proposed a novel mechanism that passive 
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fluctuations can produce PBHs with scales exiting horizon near the end of inflation. This would result in light PBHs 
with mass M < 10 15 g that may have interesting cosmological consequences in the early universe. We have discussed 
three models in which passive fluctuations can easily generate PBHs, noting that the characteristic power spectrum 
is rather generic, namely, the spectrum is extremely blue at small scales. 
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